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ABSTRACT 


A wave height transformation model is developed using surface roller theory. Roller energy production 


is included in the energy flux balance equation to predict rms wave height for randomly varying, irregular 


waves over arbitrary oathymetry. The dissipation function is defined using wave roller theory, where the 
area of the roller is defined from a simple bore analogy. The Rayleigh distribution is used to statistically 
describe wave heights as waves shoal, break, and dissipate. Model predictions are compared with data 
acquired on both barred and near planar beaches. The surface roller wave height transformation model 
predicts rms wave heights with an average mms error of 6.5% for a barred beach over ihree days, 3.0% 


for two planar beaches over four days, and within 4.5% average error for all locations over seven days. 


The model has two free parameters, o representing the type of breaker, and y a measure of breaking wave 
saturation, also a function of beach slope. Optimal values of both parameters are chosen by model fitting. 
The model is sensitive to y, but not o. The surface roller model improves the bore dissipation model 
[Thornton and Guza, 1983] by decreasing the average rms error by 40% while decreasing model 


sensitivity to input parameters. 
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I. INTRODUCTION 


Wave heights gradually increase as waves shoal from deeper to shallower water. As waves propagate 
into very shallow water, they become unstable and break, generating turbulence at the surface boundary 
layer. The turbulence due to wave breaking is the primary dissipative mechanism causing wave energy 
decay in the surf zone. Dissipation rate and turbulent penetration depth are dictated by breaker type, 
commonly described as spilling, plunging, or surging. While plunging wave turbulence can penetrate to 
the bottom, spilling breaker turbulence is primarily confined to the surface layer between the crest and 
trough. Two important parameters which distinguish breaker type are wave steepness and beach slope. 

Cacina, [1989] concluded that wave breaking is a function of both deep water wave steepness (H,/L,) 
and beach slope (tan), and is correlated with the deep water surf similarity parameter described by 


Battjes [1974] as 


(1) 


where H,,,,. is the deep water wave height, and L, is the deep water wave length (gT’/2r). 
Several methods have previously been employed to describe wave height distributions in the surf zone. 
The earliest models describe shoaling waves as monochromatic and entirely dependent on local water 


depth 


H_. = yd (2) 


where y is an adjustable coefficient, and d is local water depth. More recent wave height transformation 


models use the concept of energy flux balance and a physically based dissipation function. In the 
stochastic models of Battjes and Janssen [1978] and Thornton and Guza [1983] (TG83), wave heights are 
described statistically. Energy dissipation due to shallow water wave-breaking is modeled with simple 
periodic bores. 

The objective of this paper is threefold: to improve the previously developed bore dissipation model 
by Thornton and Guza, {1983] by incorporating surface roller theory; to determine the influence of wave 
rollers on wave height transformation in the surf zone; and to reduce sensitivity to quasi-arbitrary model 
parameters. Waves are modeled as random, having irregular amplitude, and approaching from angles over 
arbitrary bottom profiles. Roller mode! generated wave heights are compared with observed wave heights 
acquired on a naturally barred beach during DELILAH (Duck, N.C., 1990), and on near-planar beaches 
from NSTS (Leadbetter Beach, Santa Barbara, 1978 and Torrey Pines Beach, San Diego, California, 


1980). 


Il. SURFACE ROLLER MODEL 

As the wave shoals, its height increases and the wave eventually breaks generating turbulent kinetic 
energy (tke) in a volume of water which becomes detached (separated) from the wave form and pitched 
forward down the wave face, defined as the wave roller. The surface roller is believed to play a vital role 
in wave height transformation within the surf zone [Fredsoe and Deigaard, 1992]. Svendsen and Madsen 
[1984] first modeled breaking waves as rollers, i: which they envisioned a mass of water pushed by the 
wave front with horizontal velocity equal to the wave celerity C (Fig 1). The shear layer between the 
water in the roller and wave motion results in the dissipation of roller tke Jocally within the wave/roller 
system. Because the surface roller propagates with the wave, roller turbulence can be the main source of 
turbulence in the inner surf zone region [Basco, 1983]. 


The energy flux balance equation is used to determine wave energy gradients across the surf zone 
2 (EC cosa) = D (3) 
af 


where E is energy, C, is the wave group velocity, & is the mean incident wave angle calculated using 
Snell's Law, and D is dissipation. As in TG83, the model assumes stationary wave conditions, straight and 
parallel bottom contours, and random waves which are narrow-banded in both frequency and direction 
(i.e., all waves are from the same direction and of a single frequency). Energy flux is contained in two 


terms, representing contributions from the wave, £,, and the roller, £,, 


E-E,+E, 


The dissipation term in Equation (1) is derived from the wave roller. Thus Equation (3) becomes 


a dae 20 oe 
By FeO 4088) + Ay Ei Comse) = €, (5) 


where €, represents roller dissipation. 
Thornton and Guza {1983} show that wave height probability density functions both inside and outside 


the surf zone are well described by the Rayleigh distribution, p(H), where 


pt = Pex * yy 6) 
Haws rms 


The ensemble averaged wave energy, <£,>, is calculated by integrating £,, through the Rayleigh 


distribution using linear wave theory 


<E,> = 


| rms 


pg [ HHH) = 78H mg (7) 
0 


where p is seawater density, z is gravitational acceleration, H is wave height, and H,,,, is the root-mean- 


squared wave height. Group velocity, C, is described by linear wave theory 


where C is phase velocity, d is depth, and & is the wavenumber associated with peak frequency fe of the 


spectrum. 
Surface roller energy is given by [Svendsen, 1984] 


2 
E = BAC (9) 
2L 


where 4 is area of the roller and Z is the wavelength of the wave. Engelund [1981] uses an analogy with 


a hydraulic jump to analytically determine roller area, 


3 
H (10) 


where o is the angle of the roller/wave stress vector, Tt, with respect to horizontal (Fig 1) . The value of 
© is a function of breaker type and is believed to vary from a maximum value when breaking is initiated 
to a minimum value (estimated to be 10 degrees) in the inner surf zone [Engelund, 1981]. Thus Equation 


(9) becomes 


p= Of? (11) 


The roller dissipation term, €,, represents work done on the roller per unit area per unit time, and is 


given by 


= 3 
€, x xe = pgfAtano = gr (12) 


The work done by the roller is accomplished by the shear stress between the rollers bottom boundary and 
the wave's surface boundary, beginning at the toe (Fig 1). Turbulence generated at the toe moves 
horizon.a'.y through the top of the wave and fills the crest area as the roller slides down the wave face 
[Bradshaw and Ferris, 1965}. Turbulent kinetic energy is removed from the roller at the same rate. 
Because the generation of turbulence and energy dissipation only apply to breaking waves, it is 
necessary to identify which waves are breaking. The average rate of energy dissipation, <e,>, (as well as 
the average roller energy flux, E,) is calculated by multiplying the dissipation for a single broken wave 
of height H by the probability of wave breaking at each height, p,(H), and integrating over all H 
{Thornton and Guza, 1983}. Breaking waves are identified by weighting the Rayleigh function for all 


waves, breaking and non-breaking 


PAH) = WAP) (13) 


where W(H) is a weighting function, and the area under the distribution p,(71) is equal to the percent of 
breaking waves. The weighting function is based empirically on acquired data which identify breaking 


waves, including both incipient breakers and rollers, and is given by 


H 
WH) = (—=y[ 1 - exp - (fy) (14) 
yd yd 


W(H)< 1 unless conditions are totally saturated after which all waves break and W(H)=1. 
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Integrating H’ in Equations (11!) and (12) through the breaking wave height distribution yields 


Cd H? 1 
<H,> = [ H°pHdH) - Se Seeyy ty 
4 H 2 (15) 
0 (ydP + ( =) 


Using Equations (15), (11) and (7) in Equation (5) yields the ensemb!e averaged energy flux balance 


a1 <2) , 0 pC *fcosa His 7 1 : 
Bx BPC mal * (So hano GaP a! eee - fir 
“yd 
3Vxp —l : ) (16) 
16d (yd) Has,> 
1+ (—>) 
yd 


used to describe wave height transformation including a surface roller. For saturation conditions W(H)=1 


(when all waves are breaking), the ensemble averaged breaker height reduces to 
<H?> = ws H (17) 


simplifying the energy flux balance equation to 


8 3¥%pCfeosa ~ VRS 53 
ar 32dtano Has) 16d Hag 


d/1 = 352 
By g PBC soe Ha + 


For saturation conditions, we must determine when the probability of all waves breaking is equal to one, 


or 


J paanaan =f wanpanaey = 1 
0 0 


Substituting Equations (6) and (14) into Equation (19) and solving for the shallow water depth limiting 


condition, saturation conditions exist when 


in 
—=— = 1.27 
yd 


and Equation (18) is used instead of (16). A simple forward stepping numerical scheme is sufficiently 


accurate [Thornton and Guza, 1983] to solve Equations (16) and (18), written 


(EC). sg (E,C)), = <e,>,Ax - EC): . (E,C,), 


Due to the equation’s nonlinear nature, (H,,,), (shoreward) must be solved iteratively given (#/,,,,),, mean 


incident angle &, and mean frequency f. 


Hil. FIELD DATA 


The surface roller model is verified by comparing data acquired during three field experiments 
conducted on two near-planar beaches of different slope (NSTS) and a barred beach (DELILAH). The 
experiments used extensive alongshore and cross-shore arrays of current meters and wave sensors. Only 
the cross-shore arrays are considered here. Days for data comparison are chosen from each experiment 
based on model assumptions of straight and parallel contours, steady state conditions, and narrow- 
bandedness in frequency and direction. Table 1 summarizes offshore wave and beach conditions for the 
days analyzed. 

A. TORREY PINES BEACH, CALIFORNIA 
The first NSTS experiment was held in November, 1978 at Torrey Pines Beach, located north of San 


Diego, California. Beach slope was approximately linear, varying between 1:50 in the surf zone to 1:20 


in the swash region. Beach face topographies were generally concave-up and beach topography was 


relatively uniform in the alongshore direction. 

A cross-shore array of 11 current meters, 4 pressure transducers, and 4 wave staffs were used to 
measure rms wave heights across the surf zone. Data were telemetered from the instruments at Torrey 
Pines to the Shore Processes Laboratory, Scripps Institution of Oceanography and digitally recorded at 
a rate of 64 samples per second. The data was then low-pass filtered and reduced to 2 samples per second 
(Thornton and Guza, 1983}. 

A wide variety of wave and weather conditions were encountered during the experiment. Significant 
offshore wave heights, H,, varied between .6 and 1.6 meters. The average peak frequency of the incident 
wave spectra, f,, was fairly constant at 0.07 Hz (T,=14s). Incident wave angles in 10-m depth were limited 
to less than 15 degrees due to offshore island shadowing and refraction [Pawka et al., 1976]. Two days, 


4 November and 10 November, are chosen for data comparison. 


B. SANTA BARBARA, CALIFORNIA 

The second NSTS experiment was conducted at Leadbetter Beach in Santa Barbara, California over 
a one month period in the winter of 1980. Bottom contours inside 6m depth were relatively straight and 
parallel. Mean nearshore beach slope varied between 1:33 and 1:16, depending on tide level. 

A cross-shore array of 16 current meters and 6 pressure sensors measured the cross-shore wave 
transformation properties. Two wave slope arrays in deep water measured incident wave spectra. All 
instrument signals were cabled to shore and digitally recorded [Thornton and Guza, 1986]. 

Leadbetter Beach is oriented approximately east-west on a predominately north-south California coast. 
To reach Santa Barbara, the open-ocean North Pacific swell must pass between Point Conception and the 
Channel Islands, a narrow window of +9 degrees centered on 249 degrees true [Thornton and Guza, 
1986). This unique beach orientation provided good data sets 3-6 February, of which days 4 and 5 are 
chosen for model comparison. 

C. DELILAH 

The third data set used was acquired on a barred beach during DELILAH held October 1-21, 1990 
at the US Army Comps of Engineer's Field Research Facility at Duck, North Carolina. Both bar location 
and beach slope, which varied between 1:25 and 1:19, responded to a wide variety of wave conditions 
[Church and Thornton, 1993}. 

Directional wave spectra were obtained in 8m depth from an alongshore array of bottom mounted 
pressure sensors. A cross-shore array of 9 pressure sensors were located from 4m depth to the shoreline. 
The array provided near-continuous data throughout the experiment with a sampling rate of 8Hz. An 
autonomous Coastal! Research Amphibious Buggy (CRAB) was used for daily bathymetric measurements. 

Three days (Oct 10, 12, and 14) are cnosen for model comparison when waves were relatively narrow- 
banded, varying in rms wave height from .77 to 1.30 meters, and arriving at angles between | and 12 


degrees from normal to the beach. 
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IV. MODEL RESULTS 


Model results for representative days during each experiment are presented in Figures 2-8. The two 
free parameters, Y and o, of the roller model are determined by fitting model H,,, to observed data in a 
least square sense. Optimal values for each day and the resultant percent mms error are given in Table 1. 
Root mean square error is determined at the wave height observation locations. The predicted and 
observed wave heights, and dissipation as a function of offshore distance and depth are shown in Figures 
2a-8a. The gradient wave and gradient roller energy fluxes as a function of offshore distance and depth 
are shown in Figures 2b-8b. Dissipation (modeled by Equation (12)) in Figures 2a-8a is the sum of the 
gradient energy fluxes in Figures 2b-8b. The percent of breaking waves as a function of offshore distance 
and depth are shown in Figures 2c-8c. 

The DELILAH results indicate two significant wave-breaking regions. The first region just seaward 
of the offshore bar is relatively broad, while the second region shoreward is more abrupt. Weather and 
wave conditions discussed earlier are evident. On October 10, relatively small waves with an average 
period of 9.1 seconds arrived at an average incident angle of 1.9 degrees. The tide level is .7m. As a 
result, only 55% of the waves break at the offshore bar. Some waves continue breaking as they propagate 
shoreward, until at approximately 20 meters offshore the remaining waves break and saturation occurs. 

Waves began building on the 1 Ith due to a remote storm and continue building through the 13th. By 
12 October, larger waves are approaching with an average period of 15.6 seconds at incident angles up 
to 12 degrees. The tide level at observation time is .8 meters. The majority of waves (80%) are predicted 
breaking at the offshore bar due to the wave height increase, resulting in additional dissipation. Some 
waves continue to break shoreward of the bar until reaching saturation at approximately 25 meters offshore 
where all waves break. Note how the bathymetric trough has widened by approximately 20 meters since 
10 October, and beach slope has decreased. By October 14, smaller waves are arriving at an average 


period of 10 seconds and average incident angle of 7 degrees. At a low tide level of -.2 meters at 


il 


observation time, saturation occurs seaward of the bar. The dissipation region has narrowed and fewer 
waves reform in the trough, now widened by the storm to 50 meters. Beach slope has further decreased 
shoreward of the bar. 

Using data acquired at 5 beaches, Sallenger and Holman, [1984] applied linear regression to conclude 


that for saturated breaking waves, y is related to beach slope, f, (Fig 9) by the relation 


Yoms = 3-2tanB + 0.30 


and that y does not appear to depend on wave steepness. Optimal and o values reflect wave and beach 
conditions (Table 1) and confirm the linear relationship in Equation (22). How these free parameters 
reflect changing conditions can be best observed during DELILAH. On October 10th, beach slope is 
steepest requiring a relatively large y. On 12 October, larger waves have decreased the beach slope, 
reflected in the y decrease. Finally, by 14 October the beach slope has decreased significantly resulting 
in a major decrease in y. 

Roller steepness, sino, also changes with wave conditions; the steepest angles are required for 
plunging breakers, such as those occurring on the 12th, and smaller angles for spilling breakers. The same 
trends hold for the near-planar Leadbetter and Torrey Pines beaches, but with smaller percent nns error 


due to constant beach slopes. 


V. DISCUSSION AND CONCLUSIONS 


Sensitivity tests varying o and y versus rms wave height prediction error are shown in Figures 10a 
and 10b. All model runs, regardless of wave or beach conditions, produce the same trends. Roller 
steepness can vary up to 40 degrees in most cases negligibly changing rms error. The reason for 
insensitivity to 6 is seen by examining the energy flux balance Equation (16). 

Roller steepness is used to describe both the surface roller production and dissipation terms. After 
substituting tano into Equation (10) into the dissipation term Equation (12), tano cancels and remains 
only in the roller energy term. As shown in Figures 2b-8b, the gradient wave energy flux is more than 
two orders of magnitude larger than the gradient roller energy flux term. This effectively reduces the 
surface roller model free parameters from two to one in y. 

The model is sensitive to changes in gamma. The expected range of y using Equation (22) is .3 to 
.6 for the range of beach slopes from 0-.1. Examples of varying gamma for planar and barred beaches are 
shown in Figure 1} 1b. Model error is fairly insensitive for the barred beach (42%) when gamma is within 
the range .3 to .S. Error rapidly increases for gamma outside this range. The reason for the rate of change 
difference can be seen analyzing small and large gammas in Equation (16). As rms wave height is raised 
to the Sth power in the roller energy and dissipation terms compared with the 2nd power as in the wave 
energy term, a small decrease in causes the magnitude of the gradient roller energy flux and dissipation 
terms to rapidly approach wave energy flux magnitude. This rapid change in magnitude causes rapid 
changes in rms wave height predictio’: and percent rms error. A large positive change in Y only decreases 
the roller energy and dissipation terms more, which are already small relative to the wave term, resulting 
in slower percent mms error rate changes. 

Relatively low values of y compared to the expected range (.3-.6) are required to minimize the rms 


error for Torrey Pines data (Table 1). The reason is hypothesized to be linear bore theory deficiencies 
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based on conclusions by Cacina, [1989]. A review of the bore dissipation model is required for further 


investigation. 

The bore dissipation model [Thornton and Guza, 1983] predicts shoreward rms wave heights within 
an rms error of +9% for random waves and arbitrary bottom profiles using the same data for Torrey Pines 
and Leadbetter Beaches. The bore model energy flux balance equation is given by Equation (3), D being 
described by linear bore dissipation. Frictional dissipation at the bottom boundary layer is shown to be 
negligible (less than 3% of breaking wave dissipation within the surf zone) except in the run-up region 
and it is subsequently neglected. 

Applying conservation of mass and momentum at regions of flow upstream and downstream of the 


bore, the average rate of energy dissipation per unit area is calculated [Stoker, 1957] 


Oe A Ms 
S408 Q = 798 7. Q (23) 


where B is a breaker coefficient of O(1), representing the fraction of foam on the wave face (accounts for 
various breaker types). Q is the volume discharge per unit area across the bore (Fig 11), described simplest 


by Hwang and Divoky, [1970] as 
Ch 
2 am 24 
Q L (24) 


where C is wave speed and L wavelength. Substitution into Equation (23) yields 
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The bore dissipation function has the same form as roller dissipation with the exception of the B 


coefficient. The reason for their similarity is because a hydraulic jump (or bore) is used to calculate the 
area of a roller (Equation 10) when deriving roller dissipation. The fact that they both are similar after 
derivation from two perspectives would lead one to expect similar results. 

The surface roller model differs from the bore dissipation model by only the roller energy term. 
Although the roller term is much smaller than the wave energy term, it is a significant addition to the bore 
model. The term provides two changes. First, it allows the propagation of dissipation shoreward. Bore 
model dissipation results for 12 October (DELILAH) and 4 February at Santa Barbara (Fig 12) are 
compared with the roller model (Figs 3a and Sa). The dissipation plots are similar for both models at the 
planar beach, and in the offshore bar region at Duck. The significant difference occurs in the trough at 
the barred beach. In the bore model, waves stop breaking in the trough moving dissipation towards zero. 
In the roller model, waves break at the offshore bar and the roller continues through the trough and 
nearshore region, which is more in line with observations. The net result is the roller model more 
accurately predicts shoreward wave heights while providing more accurate offshore dissipation values, as 
shown by the numerical scheme in Equation (21). 

Secondly, the roller model decreases the number of sensitive input parameters. Thornion and Guza, 
[1983] conclude simple bore theory underestimates dissipation, for which the B parameter compensates. 
They also conclude that B and y could be combined into one coefficient, but are left separate for greater 


physical insight. Cacina, [1989] combined the B and ¥ term, called the BG parameter, described by 


1S 


The bore dissipation mode] was iterated over the BG parameter to reduce percent mms error in a least 


squares sense. Results of one barred beach day are shown in Figure 13, revealing a similar free 
parameter/rms error curve as the roller model. 

Cacina, [1989] concluded that B, a measure of wave breaking intensity, is correlated with the deep 
water surf similarity parameter described by Equation (1). It should follow that the roller model's single 
sensitive free parameter, y, be a function of beach slope and the deep water surf similarity parameter. 
Results of varying with these parameters for all days over various bottom profiles are shown in Figures 
14-16 and Table 1. An approximate linear relationship exists between 7, tanB, and ¢, confirming previous 
conclusions by Battjes {1974}, and Sallenger and Holman [1984). A high correlation coefficient of 0.93 
is found between y and tanB (Fig 14). Gamma is not as highly correlated (.67) with wave steepness 
(H,/L,) (Fig 15) or the deep water surf similarity parameter (.58) (Fig 16) . 

Roller model accuracy may be improved by using variable coefficients across the surf zone. Constant 
Ys and o's were used for each run. Optimal values were chosen by fitting model generated H.,,, values 
to observed wave heights. Improvement would require y and o be specified throughout the wave 
prediction region. 

In applying the model based on this analysis, using y=.31 gives wave height prediction accuracy 
within 10% rms error for all beach profiles. For barred beaches, a value of y=.34 is suggested; for planar 
beaches, use 7=.30 (averages for barred and planar beach days). o=20 degrees is sufficiently accurate for 
all profiles. The lower than expected values of y according to Equation (22) indicates that the roller model 


¥ is accounting for other factors associated with complex wave breaking processes. 
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The roller model improves the bore model by decreasing average mms error from 8.6% to 4.5% for 


all breaker types and variable bathymetry, accomplished by more accurately describing dissipation 
throughout the surf zone. The ruuer model simultaneously reduces the requirement of two free parameters 


to essentially one in ¥. 
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TABLE 1. ROLLER MODEL PARAMETER RESULTS. Free parameter values are given with 
resulting predicted vs observed wave height rms error. 1. Beach slope- mean slopes for planar beach days; 
foreshore slopes for barred beach days. 2. Mean incident wave frequency in Hertz. 3. Deep water wave 
steepness parameter.4. Deep water surf similarity parameter. 5. Mean incident wave angie in degrees. 6. 
The number of wave height observations for a given day. 7. Roller steepness angle- an adjustable free 
parameter in the roller model. 8. Gamma- an adjustable free parameter in the roller model. 9. Percent root 
mean squared error for observed vs model predicted wave heights. 


pats tenp} = ¢ (Hs)? H,/L, ([.° as #obs* o’ y° =‘ t8rror’ 
100ct 053 103 -0055 719 1.9 07 30 37 4.7 
4 

120ct 042 064 .0027 808 11.3 06 40 34 5.1 
140ct 040 103 .0049 582 6.5 08 30 30 9.7 
04Feb 038 070 0024 875 9.0 13 5 32 2.6 
OSFeb 035 078 .0023 801 8.4 10 20 32 2.5 
04Nov 021 070 0010 662 0.0 12 10 28 3.7 
10Nov 021 063 0014 §72 0.0 11 1s 27 3.1 

a NN Ta a a 
All days 67 21 31 4.5 


C = phase speed of wave 
A = area of roller 

L = length of jump 

6 = roller steepness angie 
h = mean water depth 

H = height of wave/bore 


Figure 1. Surface Roller on a wave face depicting stress vector (t) and roller steepness angle (¢) 
with respect to horizontal. 
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Figure 2a. Roller model wave height prediction using y=.37, 0=30 degrees. 
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Figure 2b. Roller model gradient wave and gradient roller energy flux. y=.37, 
o=30 degrees. 
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Figure 2c. Percent of waves breaking as a function of offshore distance and depth. 
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Figure 3a. Roller model wave height prediction using y=.34, o=40 degrees. 
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Figure 3b. Roller model gradient wave and gradient roller energy flux. y=.34, 
o=40 degrees. 
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_ Figure 3c. Percent of waves breaking as a function of offshore distance and depth. 
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Figure 4a. Roller model wave height prediction using y=.30, 0=30 degrees. 
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Figure 4b. Roller model gradient wave and gradient roller energy flux. ¥=.30, 


o=30 degrees. 
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Figure 4c. Percent of waves breaking as a function of offshore distance and depth. 
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Figure Sa. Roller model wave height prediction using y=.32, o=5 degrees. 
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Figure 5b. Roller model gradient wave and gradient roller energy flux. y=.32, 
o=5 degrees. 
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Figure Sc. Percent of waves breaking as a function of offshore distance and depth. 
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Figure 6a. Roller model wave beight prediction using y=.32, o=20 degrees. 
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Figure 6b. Roller model gradient wave and gradient roller energy flux. y=.32, 
o=20 degrees. 
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Figure 6c. Percent of waves breaking as a function of offshore distance and depth. 
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Figure 7a. Roller mode! wave height prediction using y=.28, o=10 degrees. 
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Figure 7b. Roller model gradient wave and gradient roller energy flux. ¥=.28, 
o=10 degrees. 
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Figure 7c. Percent of waves breaking as a function of offshore distance and depth. 
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Figure 8a. Roller model wave height prediction using y=.27, o=15 degrees. 
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Figure 8b. Roller model gradient wave and gradient roller energy flux. y=.27, 
o=15 degrees. 
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Figure 8c. Percent of waves breaking as a function of offshore distance and depth. 
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Figure 9. Gamma and beach slope relationship (Sallenger and Holman, 1984). 
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Figure 10a. Gamma sensitivity test using optimum o's from Table 1. 
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Figure 10b. Sigma sensitivity test using optimum ys from Table 1. 
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Figure 11. Periodic bore used to describe spilling breakers (from Thornton and 
Guza, 1983). 
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Figure 12. DELILAH experiment. Bore Model rms wave height prediction in meters (solid line above zero 
reference line), turbulent kinetic energy calculation or dissipation (dashed line), zero reference (dotted 
line), and depth in meters (solid line below zero reference) on days indicated (Church and Thornton, 
1993). 
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Figure 13. Variation in rms error due to change in BG parameter for a bore dissipation model run on a 
barred beach during SUPERDUCK experiment in Duck, N. Carolina (from Cacina, [1989)). 
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Figure 14. Correlation of y and beach slope, tanp. 
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Figure 15. Correlation of y and deep water wave steepness (H,/L,). 
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Figure 16. Correlation of y and the deep water surf similarity parameter. 
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